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A generalized definition of intrinsic and extrinsic transport coefficients is introduced. We sfiow 
that transport coefficients from the intrinsic origin are solely determined by local electronic structure, 
and thus the intrinsic spin Hall effect is not a transport phenomenon. The intrinsic spin Hall current 
is always accompanied by an equal but opposite intrinsic orbital-angular-momentum Hall current. 
We prove that the intrinsic spin Hall effect does not induce a spin accumulation at the edge of the 
sample or near the interface. 

PACS numbers: 72.10.-d, 72.15.Gd, 73.50.Jt 



Recently, there are emerging theoretical interests on 
the spin Hall effect in a s pin- orbit coupled system 
laaiiaiillSiliailllllllllll The splnHaU 
effect refers to a non-zero spin current in the direction 
transverse to the direction of the applied electric field. 
Earher studies had been focused on an extrinsic effect 
[THl ITfil |. namely, when conduction electrons scatter off 
an impurity with the spin-orbit interaction, the electrons 
tend to deflect to the left (right) more than to the right 
(left) for a given spin orientation of the electrons. Thus 
the impurity is the prerequisite in the extrinsic spin Hall 
effect. Recently, the spin Hall effect has been extended to 
semiconductor heterostructures where the spin-orbit cou- 
pled bands are important. It has been shown that the spin 
current exists in the absence of impurities, termed as the 
intrinsic or dissipationless spin Hall effect (ISHE) in or- 
der to distinguish the impurity-driven extrinsic spin Hall 
effect (ESHE) mentioned above. In general, the magni- 
tude of ISHE is two to three orders larger than that of 
ESHE; this immediately generates an explosive interest 
in theoretical research on the ISHE since the spin cur- 
rent is regarded as one of the key variables in spintronics 
application. 

However, the spin current generated via ISHE is fun- 
damentally different from conventional spin-polarized 
transport in many ways. First, the spin current is car- 
ried by the entire spin-orbit coupled Fermi sea, not just 
electrons or holes at the Fermi level Hi]. Second, ISHE 
exists even for an equilibrium system (without external 
electric fields) ^ and ISHE is closely related to the dielec- 
tric response function that characterizes the electronic 
deformation jfl]. Most recently, it is proposed that the 
intrinsic spin Hall effect exists even in insulators [T7|. 
The above unconventional properties cast serious doubts 
on experimental relevance of the intrinsic spin current. 
It has been already alerted by Rashba that the 

ISHE may not be a transport phenomenon. Due to the 
ill-defined nature of the spin current in the spin-orbit cou- 
pled Hamiltonian, theories utilizing different approaches 
produce contradicting results: some predicted a zero spin 
Hall current in the presence of an arbitrary weak disor- 
der and some claimed a universal spin conductivity at 
weak disorder. In this letter, we do not try to resolve 
the above theoretical debate, instead we reveal the spu- 



rious nature of the intrinsic spin Hall effect and discuss 
its experimental consequences. We first define general- 
ized intrinsic and extrinsic transport coefficients from the 
semiclassical transport equation. We show that the in- 
trinsic spin current is always accompanied by an equal 
but opposite orbital-angular-momentum (0AM) current 
for a spin orbit coupled system. Thus, the intrinsic mag- 
netization current which is the sum of the spin current 
and the orbital angular momentum current is identically 
zero. Next, we construct the equation of motion for the 
spin density in the presence of the intrinsic and extrin- 
sic mechanisms. We find that the intrinsic spin current 
is exactly canceled by a spin torque and thus the spin 
accumulation at the edge of the sample is solely deter- 
mined by the extrinsic spin current. The above results 
make us conclude that the intrinsic spin current has no 
experimental consequences in terms of the spin transport 
measurement for an arbitrary strength of the intrinsic 
Hall conductivity. Therefore, the intrinsic spin current 
is a pure theoretical object, at least, in the limit of the 
semiclassical picture of the spin transport. Finally, we 
brief comment on the most recent experimental results 

Let us consider a spin-dependent Hamiltonian 

H = -^y^ + V{v,a)+eEx + V,{T,(T) (1) 

where the second term represents a periodic spin-orbit 
potential, the third term is the interaction with a DC 
electric field E in the x-direction, and the last term is the 
impurity potential that may or may not depend on spin. 
Now let us consider how an arbitrary dynamic variable 
G responds to the electric field. A standard semiclassical 
version of the linear response function J^{v) is 



^(r)=^GkA/(ekA,r) 



(2) 



kA 



where G can be any dynamic variable such as the current 
density, the spin current density, the magnetic moment, 
etc., GkA = /dr«'+^(r)G*kA(r) =< *kA|G|*kA > is 
the expectation value for the eigenstate 'I'kA(r) (Bloch 
states) determined by the first three terms in Eq. (1), 
A = ±1 represents the index of the spin sub-band, and 
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/(eicA, r) is the distribution function that depends on the 
detail of the scattering potential Vi{r,a), the last term 
of Eq. (1). The dependence of J7(r) on the electric field 
enters in two places; the wavefunctions and the distribu- 
tion function. We may expand them up to the first order 
in the electric field. The wavcfunction is written as, 



(1)/ 



(3) 



where ^1^^"^^ (r) is the unperturbed electronic structure de- 
termined by the first two terms in Eq. (1), and 

-*L°l(r), (4) 



E 

k'A^kA 



k'A' 

fkA — Ck'A' 



is the first order perturbation to the third term in Eq. (1). 
Similarly, we write the distribution function in terms of 
the equilibrium and non-equilibrium parts. 



/(ekA,r) = /"(ekA) + 



(5) 



where /"^ is the equilibrium distribution function and the 
non-equilibrium function (7(kA,r) is proportional to the 
electric field. By placing Eqs. (3) and (5) into Eq. (2) 
and keeping only the first order term in the electric field, 
we have J{r) = + Jext where 



J,„, = 2Re^ < *L°a|G|vI/W > f{e^>.) 



(6) 



kA 



is defined as the intrinsic linear response and Re stands 
for the real part, and 

Je.^ = E < > {-^) 5(kA,r) (7) 



kA 



kA 



is called the extrinsic linear response. The above dis- 
tinction between intrinsic and extrinsic contributions to 
the transport properties has been already introduced by 
a number of groups, in particular, by Jungwirth et al. 
[l^ in their study of the anomalous Hall effect in itiner- 
ant ferromagnets. Equation (6) shows that the intrinsic 
linear response coefficient is not related to the transport 
phenomenon since jTint is determined by the equilibrium 
distribution function and the local electronic structure. 
Thus, there arc no transport length scales such as the 
mean free path or spin diffusion length in jTint- The ex- 
trinsic linear response, Sfext, is a true transport quantity 
because it is directly proportional to the non-equilibrium 
distribution function that is determined by various scat- 
tering mechanisms. At low temperature, the factor of 
df^ /de\^\ limits the transport states to the Fermi level. 
Comparing Eq. (6) and Eq. (7), we realize that the intrin- 
sic effect is simple and easy to calculate while the extrin- 
sic effect is much more complicated. As long as we know 
the Bloch states the intrinsic transport coefficient 

can be straight-forwardly evaluated since is known. 



The extrinsic transport coefficient not only depends on 
the Bloch states, but also on the non-equilibrium distri- 
bution function that is usually the center of the relevant 
physics. Here, however, we should concentrate on the 
easy problem: calculation of the intrinsic transport from 
Eq. (6) by using a model Hamiltonian. 

We choose a Rashba Hamiltonian to illustrate the 
physics of the intrinsic transport properties. A similar 
calculation can also be performed for a Luttinger Hamil- 
tonian ||2fli. For the Rashba Hamiltonian, the second 
term of Eq. (1) is V{y, a) ~ {a/K)cr • (p x z) where p is 
the momentum in xy plane, cr is the Pauli matrix and 
a is the coupling constant. Before we calculate the spin 
and the 0AM Hall currents from Eq. (6), we list the 
wavefunction and the dispersion relation of the Rashba 
Hamiltonian so that one can easily follow our derivation 
at each step 
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2A V ^^K^x + iky)k 



(8) 



where A is the area of the 2-dimensional electron gas, 
and 
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EkA 



2m 



Xak 



(9) 



where fc = |k| = yk-^ + ky. By placing above two equa- 
tions into (4), we have 

XeEk,, , i-n-i , s ^^^.^ 



fr) = 



In obtaining the above result, we have used < 



,(0) 



AA'. 
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We now proceed to calculate the spin Hall current by 
taking the operator G = {l/2)[szVy + VySz] where Sz = 
{h/2)az is the z-component of the spin operator and Vy 
is the y-component of the velocity. By placing the above 
definition along with Eqs. (8) and (10), and by taking the 
distribution function a step function at zero temperature, 
we obtain the intrinsic spin Hall current from Eq. (6), 



(11) 



where we have assumed that the Fermi energy is larger 
than the spin-orbit coupling energy so that both spin sub- 
bands cross the Fermi level. Equation (11) represents 
the universal spin conductivity (c/Stt) obtained by many 
groups 

Our central question is: what is the physical meaning 
of this spin current derived from the equilibrium distribu- 
tion function? To see clearly what this spin current rep- 
resents, we recall that the spin is not a conserved quantity 
in a spin-orbit coupled system. If one is interested in the 
magnetization current or the total angular momentum 
current, one should also include the 0AM Hall current. 
The 0AM Hall current can be similarly calculated by 
introducing an operator for the 0AM Hall current 



G 



1 



[(r X p) 



Vy{r X p)^ 



(12) 
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where (r x p)^ = xpy — ypx is z-component of the 0AM. 
The same straightforward evaluation of Eq. (6) leads to 



J' 



orbit 



Stt 



(13) 



Thus the 0AM current is exactly equal and opposite to 
the spin current. This result is not surprising at all: 
the total angular momentum (z-component), spin plus 
orbital, is conserved for the Rashba Hamiltonian and 
thus we can choose the Bloch states that are simulta- 
neous cigcnstates of the total angular momentum and 
the Hamiltonian. In fact, one can directly show that 
the total angular momentum current vanishes if we use 
\sz + Lz,H] = to the Rashba Hamiltonian. 

Having discussed that the intrinsic spin current is al- 
ways accompanied with the 0AM current in a bulk spin- 
orbit coupled material, our next question is whether the 
intrinsic spin current can produce a spin accumulation at 
the edge of the sample? To answer this question, we recall 
the basic idea of the spin accumulation for the extrinsic 
spin current. When an extrinsic spin current spatially 
varies, non-equilibrium spins will be accumulated so that 
the spin-diffusion is balanced by the spin-drift current. 
Equivalently, the spin accumulation results in the chemi- 
cal potential splitting between two spin sub-bands and a 
voltage can be measured experimentally when the sample 
is attached to a ferromagnetic lead . Mathematically, 
the non-equilibrium spin-dependent chemical potential or 
spin accumulation is the average of the non- equilibrium 
distribution function |2l| . For the intrinsic spin Hall cur- 
rent, the distribution is an equilibrium distribution and 
one would expect that the concept of the spin-dependent 
chemical potential breaks down. Indeed, we show next 
that the the intrinsic spin Hall current does not lead to 
spin accumulations at the sample edge and across an in- 
terface. 

To calculate the spin accumulation or the position- 
dependent spin density S{r,t) at the edge of the sam- 
ple or across an interface, one relies on the semiclassical 
equation of motion that can be generally written as 



dS{r,t) 
dt 



where Jint and Jext are the intrinsic and extrinsic spin 
current densities, r and Text are the intrinsic and extrin- 
sic spin torques due to non-commutivity of the Hamilto- 
nian with the spin operator, i.e., the spin torque is cal- 
culated by replacing G by [s,H]/ih in Eqs. (6) and (7), 
and the last term in Eq. (14) is a collision term that is to 
relax the nonequilibrium distribution function to an equi- 
librium one. To explicitly obtain the spin accumulation 
S(r, i) in a closed form, it is necessary to use a wave- 
package description so that the position-dependence can 
be readily included. Culcer et al. have already formu- 
lated that the spin torque can be written as two terms. 
In our notation, we find the expression for the intrinsic 



spin torque is r = Tq + Ti where 

To = 5] <kA|i:[s,i7]|kA>/0(kA,r) (15) 



kA 



and 



n = -V • ^ < kA|l[s, i/]r|kA > /0(kA, r) (16) 

kA 

where the symmetrization of the product of [s, H] and r is 
implied. We emphasize that Ti comes from the position- 
dependence of the center of the wave-packet. By using 
the fact that the wavefunction is an eigenstate of the 
Hamiltonian, i/|kA >= ii^kA|kA >, we immediately see 
that the expectation value of [s, iJ] is zero, i.e., tq — 0. 
To calculate ti, we use the commuting relation r] — 
—ih\'+ia{ez x cr), where v is the velocity operator. After 
a straight-forward algebra simplification, we have found 

M 



r = V.^<kA| — 

kA 



■ SV 



|kA>/°(kA,r) = V-J™t. 

(17) 

Therefore, the intrinsic spin torque exactly equals the 
divergence of the intrinsic spin current. The equation of 
motion, Eq. (14), now becomes 



dSir,t) 
dt 



V J, 



dS 

'at 



(18) 



We conclude that the intrinsic spin current does not en- 
ter into the play in the equation of motion. The spin 
accumulation S(r, t) is solely determined by the extrinsic 
part of the current density je^:*, the spin torque Text, and 
the spin relaxation in the collision term. 

We now return to our central issue on the problem, 
namely, whether the intrinsic spin Hall conductivity can 
be measured via conventional meanings of the spin trans- 
port. Since the spin current is not directly measurable, 
two schemes are usually employed: one is the realization 
of measuring the electric field induced by the magnetiza- 
tion current 24] and the other is the spin accumulation 
at the sample of the edge or across an interface [l^ ■ We 
should discuss them separately below. 

If there is a net magnetization current in a bulk ma- 
terial, a circular electric field outside the sample will be 
induced. This phenomenon is analogous to the magnetic 
field induced by a charge current, known as Biot-Savart 
law or Ampere's law. For example, it was proposed 
that a spin current generated via spin waves propaga- 
tion through a nanowire can be detected by an induced 
electric field just outside the nanowire However, 
the above proposal is applied to the case where the 0AM 
current is absent. In the present case, the 0AM current 
is exactly opposite to the spin current so that the net 
magnetization current is zero. Therefore, we conclude 
that there is no electric field associated with the intrinsic 
spin current. 
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The more efficient method to detect the spin current is 
by measuring the spin accumulation due to spatial vari- 
ation of the spin current, e.g., the Johnson-Silsbee's ex- 
periment '23'! . Based on the equation of motion given by 
Eq. (18), the divergent of the extrinsic but not intrinsic 
spin currents can lead to a buildup of spin accumulation. 
To determine the spin accumulation, one usually makes a 
relaxation-time approximation so that the collision term 
in Eq. (18) is modeled by —S/tsj. Since the intrinsic 
spin current does not contribute to the equation of mo- 
tion for the spin accumulation, the measurement based 
on the detection of the spin accumulation will produce 
a null contribution from the spin Hall effect, no matter 
how large the intrinsic spin current is. 

Two experimental groups have recently observed the 
spin Hall effect by detecting spin accumulation at the 
edges of the samples Kato et al. argued that the 

effect is extrinsic based on their experimental results that 
the spin accumulation is independent of the strain direc- 
tion, while Wunderlich et al claimed that their observed 
effect is intrinsic based on the assumption that the im- 
purity scattering is weaker than the spin-orbit coupling 
in their samples (clean limit). We point our here that 
the clean limit in the experiment does not imply the spin 
accumulation from the intrinsic origin. Instead, our anal- 
ysis has shown that no matter how large is the intrinsic 
spin current, the observed effect has to be an extrinsic 
origin because the spin accumulation is independent of 
the intrinsic spin current. 

We finally draw a picture on why the intrinsic spin 
Hall fails to produce experimental consequences. Con- 
sider a contact between a Rashba material and a non- 
magnetic material with no spin-orbit coupling. The spin 



current, as well as the orbital angular momentum cur- 
rent, would exist in the the Rashba material. However, 
the spin current drops to zero across the interface of the 
non-magnetic material, i.e., the spin current is not con- 
tinuous; this is because the spin torque produces a mech- 
anism to transfer the spin current to the orbital angular 
momentum current or vice versus. As a result, when the 
spin orbit coupling vanishes at the non-spin-orbit coupled 
material, both the spin and orbital angular momentum 
currents drop to zero. The loss of the spin current ex- 
actly equals to the gain of the 0AM current so that the 
total angular momentum current or magnetization cur- 
rent is continuous across the interface of the layers; they 
are both zero. For the same reason, the edge of the sam- 
ple never develops spin accumulation because the usual 
boundary condition of zero spin current at the surface 
is no more valid, instead, the total angular momentum 
current is zero at the surface for the intrinsic spin Hall 
effect. 

In conclusion, we have constructed a general frame- 
work for calculating intrinsic linear response coefficients. 
We have shown that the intrinsic spin Hall effect is ac- 
companied by the intrinsic orbital-angular-momentum 
Hall effect so that the magnetization current is zero in 
a spin-orbit coupled system. The intrinsic spin Hall ef- 
fect is not a useful source of spin currents because the 
intrinsic spin current does not enter into the equation of 
motion for the spin transport. Most of the proposed ex- 
perimental detections of the intrinsic spin Hall effect are 
the artifact of the boundary conditions that are not valid 
for the intrinsic spin Hall current. 

This work is partially supported by DARPA-SPINS 
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